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Malaria resurgence in African highlands in the 1990s has raised questions about the underlying drivers of 
the increase in disease incidence including the role of El-Nino-Southern Oscillation (ENSO). However, 
climatic anomalies other than the ENSO are clearly associated with malaria outbreaks in the highlands. Here 
we show that the Indian Ocean Dipole (IOD), a coupled ocean-atmosphere interaction in the Indian Ocean, 
affected highland malaria re-emergence. Using cross-wavelet coherence analysis, we found four-year long 
coherent cycles between the malaria time series and the dipole mode index (DMI) in the 1990s in three 
highland localities. Conversely, we found a less pronounced coherence between malaria and DMI in lowland 
localities. The highland/lowland contrast can be explained by the effects of mesoscale systems generated by 
Lake Victoria on its climate basin. Our results support the need to consider IOD as a driving force in the 
resurgence of malaria in the East African highlands. 

Malaria resurgence in the East African highlands in the 1990s caused serious mortality and morbidity 1,2 . 
Although climate change in general and the El-Niho-Soufhern Oscillation (ENSO) in particular were 
strongly suspected to be related to epidemics in several highland regions in early 1998 3,4 , other large scale 
climatic factors are clearly associated with malaria risk in the highlands, because epidemics also occur during 
weak ENSO events. The Indian Ocean Dipole (IOD) causes above normal rainfall in East Africa during positive 
IOD events 5,6 . The availability of mosquito breeding grounds and vector populations both increase with rainfall 7 . 
A short-term temporal association between IOD and malaria incidence in the western Kenyan highlands has been 
reported 8 . However, it is still not clear if elements of the local climate relevant for malaria transmission, specially 
rainfall 9 and malaria itself, have non-stationary signatures of the IOD on their dynamics over different time scales. 
This is of special interest provided the plausibility of regime shifts in climate caused by changes in the IOD 10 . 

East African highland environments can be expected to be more sensitive to climatic variability because they 
are farther away from water basins such as lakes, and proximity to large water bodies generally stabilizes the 
climate. However, whether the pattern of exacerbated malaria transmission observed in the Western Kenyan 
highlands and associated with the IOD, is also found in adjacent lowland populations in the Lake Victoria basin 
remains an open question. Although it has been argued that the low immunity of highland human populations 
can explain the increased size of the epidemics in the late 1990s 11 , the sensitivity of parasites and vectors to 
climatic variability can also affect the ecology of malaria transmission and more stable environments are less 
prone to epidemics 12 . Here we examine the hypotheses that 1) the re-emergence of highland malaria in the 1990s 
is associated with IOD activity and 2) the increased sensitivity to climatic variability in highland malaria trans- 
mission is underpinned by differences in environmental stability related to proximity of the various regions to 
Lake Victoria. 

Results 

Time-series of malaria hospital admissions show that epidemics began in the early 1990s and peaked in the mid to 
late 90s in the highlands (Kericho, Kisiiand Kapsabet), while epidemics began in early the 1980s and peaked in the 
mid 80s in the lowlands (Maseno and Kendu Bay) (Figure 1). In general, the association with dipole mode index 
(DMI) is stronger and wider for all the studied time scales compared with Nino3, the ENSO index (Figure 2). 
Four-year long coherent cycles between malaria and DMI were seen in the 1990s data from the three hospitals in 
the highlands (Kisii, Kapsabet and Kericho). Coherence with DMI was less obvious in the lowlands (Maseno and 
Kendu Bay) while coherence with Nino3 was more obvious in Maseno and Kendu Bay in the data from the 1980- 
90s (4-16-year cycle). These patterns are supported by cross-correlation functions (Figure SI) and by autore- 
gressive models that consider both DMI and Nino3 as forcing factors (Table SI). The best models for Kisii and 
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Figure 1 Malaria time series. Clinical records of malaria infections (A) Maseno (May 1935, November 2009, 0 00'S, 34°36'E, Altitude = 1500 m),blue 
indicates imputed values; (B) Kendu Bay (January 1980, November 2006, 0°24'S, 34°39'E, Altitude = 1240 m); (C) Kisii (January 1986, December 2000, 
0°40'S, 34°46'E, Altitude = 1670 m); (D) Kapsabet (January 1980, December 1999, 0°12'N, 35°06'E, Altitude = 2000 m); (E) Kericho (April 1965, 
November 2006, 0°23'N, 35°15'E, Altitude = 2000 m); (F) DMI, dipole mode index (dotted line) and Nino3, ENSO index (solid line, March 1958, 
December 2008). 
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Figure 2 | Cross-wavelet coherence of the malaria time series with the DMI (dipole mode index) and with Nino3 (ENSO index). (A) Maseno and DMI 
(March, 1958, December 2008), (B) Maseno andNino3 (March 1958, December 2008); (C) Kendu Bay and DMI; (D) KenduBayandNino3; (E) Kisiiand 
DMI; (F) Kisii and Nino3; (G) Kapsabet and DMI; (H) Kapsabet and Nino3; (I) Kericho and DMI; (J) Kericho and Nino3. The coherency scale is from 
zero (blue) to one (red). Red regions in the plots indicate frequencies and times for which the two series share variability. The cone of influence (within 
which results are not influenced by the edges of the data) and the significant coherent time-frequency regions (p < 0.05) are indicated by solid lines. The 
procedures and software used are those described in Chaves and Pascual (2006)™. A smoothing window of 6 months was used to compute the cross- 
wavelet coherence. 
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Kericho used a DMI with lags of 1 and 2 months, respectively, as a 
significant covariate. The local rainfall with a short-term lag (0 to 6 
months) was a significant covariate in the best models for all 5 sites. 
During the model selection process, Nino3 was discarded as a sig- 
nificant covariate for all sites. A seasonal coherence between malaria 
and rainfall was observed in all 5 sites (Figure S2) and rainfall at all 
study sites (Figure S3) had a non-stationary relationship with both 
DMI and Nino3 (Figure S4). 

Discussion 

The finding that interannual cycles of the DMI coincide with malaria 
resurgence provides strong evidence that the IOD had an effect on 
the re-emergence of highland malaria in Western Kenya, probably 
through the impact of the DMI on rainfall. This association was less 
pronounced in adjacent lowlands closer to the Lake Victoria basin. 
Mesoscale systems generated by the lake have been shown to con- 
tribute significantly to climate variability over the lake basin 13 while 
rainfall in the highlands was affected more by large-scale moisture 
from the equatorial Indian Ocean 14 . Rainfall in the highlands re- 
duces dramatically as moisture from the equatorial Indian Ocean 
diminishes (Figure 3), indicating that highland weather is more sens- 
itive to moisture from the equatorial Indian Ocean than from the 
Lake Victoria basin. The climate of the East African highlands is 
highly sensitive to perturbations during the short rainy season 15 
and the IOD may intensify and prolong the rainy season. Because 
the temperature of the highlands increases towards February, a pro- 
longed short rainy season with warmer temperatures could enhance 
the development of mosquitoes and parasites 16,17 . As a result, the IOD 
may intensify the regional malaria transmission during the short 
rainy season. This synergy effect will be more profound in the high- 
lands where temperatures are lower, and an increase in each of the 
factors (temperature and rainfall) is known to have an impact on 
mosquito populations 16,18 . 

The results from autoregressive modelling also suggest that the 
effects of the IOD vary between the highland sites. Kapsabet is in the 
highlands (altitude =2000 m) but its rainfall may be less affected by 
moisture from the equatorial Indian Ocean than the simulated rain- 
fall in Kericho and Kisii. This is reinforced by the disappearance of 
the DMI in the best autoregressive model for Kapsabet (Table SI). 
Differences between the highland sites may be due to the complex 
topography of the regions that regulates the amount of moisture 
entering from the Indian Ocean. 

The role of the ENSO in the interannual variability in East African 
rainfall has been extensively examined 19,20 . The results indicate that 
the ENSO has some influence on equatorial and coastal East African 
rainfall, with warm events being associated with higher rainfall and 
cold events with lower rainfall. However, the influence of the ENSO 
has decreased over the western Indian Ocean in recent decades, and 
the influence of the IOD on East African climates has increased 10 . 
Recent studies indicate that patterns of association between the IOD 
and rainfall are stronger, more regular and homogenous across East 
Africa 21,22 . An atmospheric circulation model suggested that Indian 
Ocean sea surface temperatures (SSTs) exert a greater influence than 
Pacific SSTs over East African rainfall 23 . Similarly, our model showed 
that the DMI had a stronger association with rainfall than Nino3 in 
the highland areas we studied. 

Although climatic phenomena are important drivers of malaria 
dynamics in the East African highlands, the IOD is not the only 
influence on the cycles of malaria incidence and non-climatic factors 
have also been implicated. Local enforcement of control measures 
potentially influences the seasonal and long-term variability of the 
incidence of malaria 12 . Chloroquine was effective in the suppression 
of malaria epidemics until the emergence of resistance 24,25 . Increased 
trends in human mobility between the highlands and endemic areas 
may also affect population vulnerability thereby modifying the mag- 
nitude of local relationships. Although these factors can modulate 



the reported associations with climatic factors, they are unlikely to 
vary in the same seasonal and interannual patterns as malaria. It 
therefore seems unlikely that these factors could produce the coher- 
ence identified in this study. However, the strength of the association 
may be greater or lesser according to such local conditions. 

The devastating epidemics had largely abated by early 2000s and in 
more recent years malaria appears diminishing in the east African 
highlands 26 . This may be the result of extensive distribution of insect- 
icide treated bed nets and the use of artemisinin-based combination 
therapies as the first-line drug 27 2 *. Further work to clarify the effect 
of IOD and other climate factors in more recent years would be of 
interest. 

Our results help in part to resolve the inconsistent associations 
reported between ENSO and malaria epidemics, and indicate that 
they merely reflect the association between the IOD and the ENSO. 
We propose that the IOD should be considered as the more appro- 
priate index to decipher the signature of climate changes in the 
resurgence of malaria in the African highlands in the 1990s. 

Methods 

Data. The primary dataset for this study was the monthly number of patients with 
malaria who were admitted to three hospitals (Kapsabet, Kisii and Kericho) in the 
western Kenyan highlands and two hospitals (Maseno and Kendu Bay) in the lowland 
areas surrounding Lake Victoria. The primary diagnosis of malaria was made by 
physicians at each hospital. All cases diagnosed in Kericho were confirmed 
microscopically, but only some cases from the other hospitals were microscopically 
confirmed. With the exception of the Maseno site, missing data were imputed by 
averaging incidence values for the month before and the month after the missing 
value. For Maseno, the monthly number of malaria cases between January 1993 and 
December 1996 was estimated by averaging the monthly values for the year before 
and the year after the missing values. The strength of the IOD was measured by the 
DMI, defined as the difference in SST anomalies between the western (10°S-10°N, 
50°-70°E) and eastern (10°S-0°, 90°-110°E) tropical Indian Ocean 5 . The DMI data 
were obtained from the Japan Agency for Marine- Earth Science and Technology 
(JAMSTEC) (www.jamstec.go.jp/frcgc/research/dl/iod/). The DMI values were 
calculated using the SST data from the HadlSST dataset. The base period for 
calculating anomalies was 1 958-2008. The obtained values were standardized to zero 
mean and unit standard deviation. The strength of the ENSO was measured by SST 
anomalies (Nino3) in the Pacific Ocean, using data from the NOAA Climate 
Prediction Center (www.cpc.ncep.noaa.gov). The base period for calculating the 
anomalies was 1971-2000. 

Statistical analysis. Nonstationary patterns of association in the time-frequency 
domain. Patterns of association in nonstationary time series can be studied using 
continuous wavelet transforms. Specifically, we used cross-wavelet coherence 
analysis to determine if peaks at a particular frequency and time in the disease 
correspond to similar peaks in a covariate. 

Time-domain time-series analysis. Based on their autocorrelation functions, we fitted 
autoregressive models to the time series and included non-linear trends obtained 
using a locally weighted smoothing (LOESS) model. The autoregressive coefficients of 
the best model for each site were used to pre- whiten the series of climatic covariates, 
thus ruling out the likelihood of a spurious correlation caused by a common cyclical 
structure on the time series. Pre-whitening is a process that rules out spurious cor- 
relations between two time series, by removing (filtering) any common structure 
between the two studied time series 30 . The residuals of the malaria time series models 
and the pre- whitened covariates were then used to estimate cross-correlation func- 
tions (CCF). Climatic covariates selected from the CCFs were included in models, 
with and without breakpoints for regime shifts. Finally, models were selected using 
the Akaike information criterion. 

Exploratory analysis. To explore the dynamics of the time series we performed an 
exploratory time series analysis and examined the autocorrelation function of the 
malaria outpatient time series 31 . An inspection of the plotted data showed that all 
malaria time series (y ( ) were likely to be first order seasonal autoregressive processes. 
Thus, we fitted the following seasonal autoregressive null model to the time series 
data: 

yi = /-' + <fri(jt-i- li) +<l>n(yt-i2- V) +<l>i<l>i2(yt-n- 1- 1 ) + *Trend + E l (1) 

For the highlands of Kisii and Kericho we had to use a second model because the time 
series for those sites were better described by a second order autoregressive process: 

yt = b>t-i — j«) + ^iyt~% — V) + xTrend + 8, (2) 

In both models, ft is the average value of the time series, trend is non-linear estimated 
using LOESS (the error was assumed to be independent and normally distributed), 
and r. t — N(Q,<j 2 ). We used this model to pre- whiten the time series of the following 
climatic covariates: rainfall, Nino3, and the DMI. Residuals of the model presented in 
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Figure 3 | Elevation (m) (A) and simulated rainfall (mm) over the Lake Victoria basin with large-scale moisture through the eastern boundary reduced 

by (B) 20% and (C) 50%. (Modified from Anyah et al. 2006 14 figure 10. ® American Meteorological Society. Reprinted with permission.) In all panels 
location color indicates the data available at each site; blue (rainfall); green (disease) and red (disease and rainfall). We used Kisumu rainfall data as a 
proxy for Kendu Bay and Maseno. Kisumu, like Kendu Bay and Maseno, is in the lowlands (A), with a similar rainfall regime (B and C) . Panel C shows box 
plots for rainfall. Rainfall variability is measured by kurtosis (K), and low values indicate a platykurtic distribution (one where conditions around the 
median are more variable). (S) indicates skewness. 
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(1) and the pre- whitened residuals of the climatic covariates were used to compute 
cross -correlation functions for the number of outpatients with each one of the 
climatic covariates. 
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